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1. Introduction 

The present work should be considered as an extension of an earlier one, [1] 
where heat-kernel asymptotics and functional determinants were computed for a 
scalar field on a bounded generalised cone (the simplest example of which is a ball). 
Here, the scalar field is replaced by more complicated bundles, specifically by spinors 
and forms. This will enable us later to discuss topological matters such as the 
analytic torsion as well as settling a few questions regarding the spectral invariants. 
However the emphasis here will be on the actual computations and the technical 
problems raised have, we think, elegant outcomes. A primary aim is to work in 
arbitrary dimensions with arbitrary p-forms. In this connection, one objective is 
to make contact with the work of Blazik et al [2] who have given general, explicit 
expressions for the first four heat-kernel coefficients for the form-valued Laplacian 
on manifolds with boundary. We anticipate going beyond their computations but, 
of course, within the confines of a restricted geometry. 

The structure of this paper is as follows. In section 2, we discuss the vari- 
ous types of eigenforms that exist on the generalised cone. For rapidity, we use 
the results of Cheeger [3] on the infinite cone, additionally imposing appropriate 
conditions at the cone boundary. 

Section 3 contains the general construction of the (function and heat-kernel 
coefficients. This closely follows our earlier development so that explanatory details 
can be kept to a minimum. The result is a formula for the p-form coefficients on the 
generalised cone in terms of those on its base for a particular, 'modified' coexact 
p-form. 

In sections 4 and 5 this apparatus is applied to the case of the ball where it 
is possible to find a compact expression for the relevant p-form (-function on the 
sphere in terms of Barnes (-functions. Another more suggestive form of this sphere 
(■-function expresses it as a sum of 0-form (-functions. The ball coefficients are 
determined as explicit functions of d and p via generalised Bernoulli functions. 

Rearranging the first few coefficients we obtain agreement with the general 
formulae of [2] evaluated on the ball. Going beyond this, we are able, in section 6, 
to find the expression for an arbitrary coefficient in terms of binomial coefficients 
and polynomials in the dimension. This is a natural generalisation of our previous 
results [1], and of Levitin's earlier [4], for the scalar case. Section 7 briefly looks at 
the important value of the (-function evaluated at zero argument. 

In section 8 we address the calculation of the functional determinants and some 
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explicit results are exhibited in an appendix. 

Related work has been performed by Elizalde et al [5] and we use their specific 
results as a check of our general forms. Their method of evaluation also works 
generally but involves an immediate reduction to a series of Hurwitz ^-functions and 
does not yield the functional form. The expressions for the functional determinants 
obtained in [5] contain integrals over \& and F functions. In an appendix we show 
how these may be reduced to derivatives of the Riemann ^-function. 

In section 9 we turn, rather more briefly, to spinors and apply our general 
method to rapidly reproduce and justify our earlier results for spectral boundary 
conditions. A curious vanishing theorem, noticed earlier, is proved in an appendix. 
Local conditions are also treated by the improved technique. 



2. Geometry and eigenforms. 

The geometry has been outlined in Cheeger [3] and already used in our earlier 
work, [1], so that it will not be described at length here. The cone manifold, M. 
has the hyperspherical polar metric ds 2 = dr 2 + r 2 dY? with dT? the metric on the 
base, A/", which is the r = 1 section of the cone. For simplicity, we shall assume in 
this paper that the base is closed, dM = 0, which is enough for the ball. 

The structure of the eigenforms of the de Rham Laplacian, A^, on the infinite 
generalised cone has been given by Cheeger [3] and, for convenience, we refer to [3] 
for the detailed action of A M on forms of the type f3 + dr A ui as well as for the 
necessary exterior calculus. (See also [6].) 

There are four basic types of eigenforms with nonzero eigenvalues, a 2 , 

aM(i) = J v( P )( ar ) ,M m 

V P r (d-l-2p)/2 > K 1 ) 



Jv(p-i)(<xr) ~~ M 



MW _ Jy(p-2)(<Xr) ~ M 

$P ~ r (d+l-2p)/2 ar/ \ a( Pp-2i 



(2) 



(3) 



(4) 



2 



where is a coexact p-eigenform of the intrinsic de Rham Laplacian, = dS+Sd, 
on the base M. We have the separation of variables relation 



Hp) = {M + ((d-l)/2-p) 2 ) 



2\l/2 



(5) 



with /u,(p) being the (coexact) eigenvalue of Ajy. We are assuming that v > 1 so 
that the 'negative' modes ~ J _ v do not arise. 

In addition, there are modes ('zero modes') whose M part is harmonic, /j, = 0, 



and 



with 



kM(E) ^g(p)( ar ) , TV 
P r (d-l-2p)/2 P 



lM(0) _ I vo(P) y 
P \ r (d+l-2p)/2 



dr A h 



AT 

■p—li 



(6) 
(7) 
(8) 



V E (P) = \(d- l)/2 -p\ = u E (d-l-p) 

u (p) = \(d+l)/2-p\ = u E (d-p). 

On M, types 1, i? and 3 are coexact and types 2, O and 4 are exact. (See the 
relations in [6].) 

For the bounded, generalised cone, conditions are to be set at r = 1. Absolute 
boundary conditions are [7,8] 



kM 
i— 3 



= o, 



r,i...j 



= 



(9) 



and have to be applied to the six types separately. 

Since <p^ , hp are pure Af forms, it is easily shown that types 1, 2, E and O 
satisfy Neumann (Robin) and types 3 and 4 Dirichlet conditions. Bessel's equation 
has to be used to derive this for type 3. 

More precisely the Robin conditions are (e.g. for type 1) 



d r ( r P-^)^j v{p) ( ar) ) 







(10) 



so that the Robin parameter is u = u a (p) = p — (d — l)/2, ([9] 3.1). The parameter 
for type 2 is u a (p — 1). 

Relative boundary conditions are obtained by dualising and, in the present 
context amount to 



(d r + d + 2-2p)(f>: 



M 

r,i...j 



= 0, ^ 



= 0. 



(11) 



3 



Now types 1, 2, E and O satisfy Dirichlet and types 3 and 4 Robin conditions. For 
type 3, Robin conditions read 

^(r (d+1)/2 - p J Kp _D(ar)) =0 (12) 

J\f 

with the Robin parameter u = u r (p) = (d + l)/2 — p = u a (d — p). The parameter 
for type 4 is u r (p — 1). Hodge ★-duality on M. thus interchanges the conditions (10) 
and (12) with 1 «-> 4 and 2 <-> 3 which differs from Cheeger. 

We denote the coexact degeneracies on Af by d(p) and remark that the exact 
degeneracies, d ex (p), and eigenvalues, n ex (p), are given by 

d ex {p) =d(p-l), n ex {p) = n{p- 1). (13) 

The structure of the eigenforms shows that the degeneracies of types 1, 2, 3 and 4 
are d(p), d(p — 1), d(p — 1) and d(p — 2) respectively. 

We now define, for later use, the modified coexact (-function on Af, the base 
(■-function, 

c j^ (s) = y- = M (u) 

2^)2. 2^ (/4(p) + ((d _ 1)/2 _ p) 2)' ^ 

(c/ [1] (3.3).) Including the 'zero modes', which is sometimes notationally conve- 
nient, one has 

which defines d(p) and where is the p-tYv Betti number of M. The summations, 
here and later, are over the mode labels which are not always explicitly displayed. 
In view of the relations (13), an 'exact' (-function is defined by 

c£» = c£i(«). (i6) 

Because the additional term ((d— l)/2 — p) 2 depends on p, Cex,p( s ) ^ s n °t the true 
exact (-function, but we expect it to be the relevant construct. 

Assume for the time being that the base Af has no boundary. Hodge duality 
on Af can be applied to yield the coexact relations 

d(d — 1 — p) = d(p), and fj,(d — 1 — p) = fx(p) (17) 

whence 

v(d- 1 -p) = u(p) (18) 
and therefore, for the coexact (-function, 

C£i- P («) = C^(«). (19) 
These relations will be verified later when Af is a sphere. 
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3. General constructions. 



The total (^-function on M. is, according to (1) - (4), (6), (7), a combination 
of exact and coexact (on M) contributions. Noting that the zero modes on Af give 
degenerate eigenvalues on M, one has 



where the af s are the eigenvalues of the p-form Laplacian on M for i-type modes. 
(20) can be written in terms of the coexact ^-function on M, (^(s), as 

C£ 4+ (s)=Ct A (s) + C£l 1 (s), (21) 

the inverse of which is 

C(«) = E(-l) p -' ? C + («)- (22) 
According to the method given in [10,9,1], and now for absolute conditions, 



+ d(p-l)lnJ Hp _ 1) (k) 

Af / 

(23) 

where the details of the contour 7 are in the above cited references. The first term 
is the Neumann (Robin) (types 1 and E) and the second term the Dirichlet (type 
3) part. For economy of writing, the degeneracy d(p) is introduced so as to take 
into account the .E-type zero modes on A/", cf (15). When p — > p — 1, the type 
1 contribution becomes a type 2, the type E a type O and the type 3 a type 4. 
For relative conditions, likewise, 



<£(«) =£ / ^r 2s k { d{p ~ 1} ln 



+ d(p) lnJ^fc) 

(24) 

It is amusing to check in detail Hodge duality on M, 

t% + (s) = 0_ p (s), (25) 
which fundamentally arises from the intertwining [8] 
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It is easily seen from (21) that (25) is equivalent to the statement that, under 
p — > d — p, C^( s ) of (23) turns into C^( s ) of (24). Firstly, this is readily seen 
to be the case from the relations (17) and (18), if the zero E modes are set aside. 
Secondly, consider the contribution coming from the zero E modes in (23). One 
has u a (p) = p — (d — l)/2 and the argument of the logarithm reads explicitly 

u aip)Ju E {p) + ^J'v E {p) w h ere I u o(p)I = v e{p)- us assume for the moment that 
u a (p) < 0, and then use the relation zJ' v {z) — uJ u (z) = —zJ u+1 (z) to write the 
above argument as —kJ UE ^ +1 . 

The factor (—k) does not contribute because it has no zero inside the contour, 7, 
and so the contribution is just that of J UE (p)+i- Now from (8), v E {p) + l = u E (d—p) 
and, taking into account Poincare duality on A/", (3^ = P^L p , the contribution of 
the p-form zero modes for absolute boundary conditions is seen to equal that of the 
(d — p)-form zero modes for relative boundary conditions, which was to be shown 
in order to complete the demonstration of Hodge duality on M. For u a {p) > use 
zJ' u {z) + vJ v [z) = zJ u _ 1 {z) to arrive at the same conclusion. 

One of our main objectives is the evaluation of the coefficients, A n / 2 , in the 

heat-kernel expansion, which is stated in generic form, 

00 

K M (r) = A n/2 r (n - D)/2 + A' M log r. (26) 

n=0 

Our preferred general computational formula, on any manifold, A4, of dimen- 
sion D, is 

A^ /2 = T((D-n)/2)ResC M ((D-n)/2), n < D. (27) 

The ability, by choosing D sufficiently large, to work with just this formula 
in order to determine any coefficient, has an additional advantage because a finite 
number of 'extra' modes does not affect the analytic structure, as in the difference 
between ( and ( when (27) is applied to the base N '. 

The calculation of the (function on M in its relation to that on A/", follows pre- 
cisely the previous pattern, [1], the Robin and Dirichlet cases now being combined. 
It leads to the following basic equation for the coexact heat-kernel coefficients on 
M. in terms of those on A/", (n < d), 



1 

+ 7«-i)/ 2 (p)-<-i)/ 2 (p-1)) 



4^(n-l)/2W - ^(n-l)MP - -UJ (28) 
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for absolute conditions. 

The A M are the heat-kernel coefficients corresponding to the base ^-function, 
(15), and the P i are known polynomials arising from the asymptotic expansion of 
Bessel functions, e.g. 

P( x )-Tx h r ((^-" + ^ 2 + & ) (29) 

" x *,b r {{D _ n + - )/2) r(6 + - /2) • l^J 

The x's correspond to Dirichlet and the z (p)'s, which depend on p through 
the u a (p), to Robin conditions. 

The combination (21) leads to the total coefficients on M, 

<1 + / 2 (P) = <i/ 2 (p) + A^ n/2 ( P - 1) (30) 

and similarly for relative conditions. (30) can be inverted as in (22). 

In general C^( s ) has a pole at s = which translates into the logr term in 
the expansion (26), the coefficient being 

A' M (p) = -Res C P M (0). (31) 

These are the only general equations that are needed but the algebra can be 
checked by confirming Hodge duality on jW, in the coefficient form 

<%(d+l-P)=<„yp)> ( 32 ) 
using the easily verified formula, 

4£(d-p) = ±4£(p) 

where, the coefficients are those resulting from the combinations of the 

coexact and 'exact' ^-functions on M discussed above, 

(t r± (s) = (i f (s)±(i r -i(s) (33) 
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4. The ball and sphere (— functions. 



We now consider forms on that cone whose base is a unit d-sphere, d = D — 1, 
i.e. on the D-ball. This will be our main application of the general formalism 
exposed in the previous sections. As we have seen, it is sufficient to look at coexact 
forms (i.e. transverse, antisymmetric fields). 

The form ^-functions (14) on the d-sphere are needed in order to find the 
corresponding heat-kernel coefficients for substitution into the fundamental equa- 
tion, (28). The spectral properties have been known for some time [11-16] and the 
coexact p-form eigenvalues of the de Rham Laplacian are readily established to be 

fi(p, /) = (/ + (d - l)/2) 2 - ((d - l)/2 - p)\ I = 1, 2, . . . . (34) 

As anticipated, we again witness the important simplification of the Bessel 
function order, (5), to the p-independent form 

u(p,l) = 1 + (d- l)/2 > (d- l)/2 (35) 

exactly as in the scalar case. The first consequence is that the absolute (^-function 
(23) simplifies to 



C» =£ / ^^^(P)^ lnJ„(*r))' 



+ d(p-l)]nJ v (k) ). (36) 

r=l 



For the time being let us work at a fixed p. The sphere coexact p-form degen- 
eracy is 

d( P i) = (2i + d-m + d-iy. 
{Pj } p \(d-p-i)\(i-i)\(i+p)(i + d- P -iy { } 

We note the symmetry, d(p, I) = d(d — 1 — p,l) and that d(d, I) = 0. In addition, 
there is a zero mode for p = and one for p = d. 
Rewrite (37) as 

(l + d-l)\ ( I 1 \ 

d ^ l > ~ p!(d_p-i)!(i_i)! \TT P ' + l + d-p-1 ) ' (38) 



and consider, firstly, the sum 



^ (l + d-l)\r l +( d -^ 
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with r = exp(— r) < 1. The idea is that this gives the square-root heat-kernel, and 
the sphere ^-function, (14), follows by Mellin transform on r as in our earlier works, 
[17,18]. 

The generating function for a given form order can be rewritten using the 
identity 

- (l + d-l)l _H_ = * (m-1)! r__ 

(1-1)1 l+p [ P ) -J^ +l (m-p-l)\(l-r)™ [ } 

which follows easily from recursion. 

There is still an overall factor of r^ d ~ 1 ^ 2 in (39) and performing the Mellin 
transform produces a series of Barnes (functions, giving, after the addition of the 
p — > d — p — 1 term, the modified coexact (function, (14), on the sphere as 



<f (*)= E ( m v 1 )c B (2 S ,(rf+i)/2|i m ) 



m=p+l 
d 



+ E („" 'Wiii,(J+l)/2|lJ 

m=d—p 



(41) 



for < p < d. Obviously (s) = 0. 

(41) is formally much simpler than expanding the degeneracy to give a series 
of Hurwitz (^-functions (a series noted in passing by Copeland and Toms, [15] and 
frequently used). This will come later. 

When p = 0, (41), with the zero mode included according to (15), gives the 
known scalar expression. The first sum reduces to a term that is cancelled by the 
zero mode and to a single Barnes (^function, the result being, [18], [1] eqn. (5.7), 

Cf = C B {2s- (d + l)/2 | l d ) + Q B {2s- (d - l)/2 | l d ). (42) 

This is elegantly seen from a rearrangement of (41), detailed in Appendix A, 
that yields (p d (s) as a finite series of scalar ^-functions by means of a recursion 
relating a p-form in d dimensions to a (p — l)-form in d — 2. The series is (0 < p < 
(d-l)/2) 

V ' (43) 



+ Vo(-7-) "("I) 



d-l\~ 2s . ^ v fd-2p-V ~ 2s 
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and this is one of our basic equations. Duality, (19), can be used to extend the 
range of p. 

When d is odd we have from (42) the special values for q e Z, 

C s Vq)=0, q>0, (44) 

which corresponds to the fact that for scalars the relevant operator is the (improved) 
geometrical Laplacian. 

From (44) follow the special values of the modified coexact function (0 < p < d) 

Cfi-Q) = Cf (-9) = (-l) P+1 ((d- l)/2 -pf\ (45) 

In particular 

Cf (0) = (-l) P+ \ (46) 

and 

C(?-i )/2 Hz)=0, q>0. (47) 

Hodge duality, as expressed in (19), can be seen explicitly from (41). Inclusion 
of the zero modes obviously violates this duality. 

It is also worth pointing out that the + combination of coexact and exact (- 
functions in (33) on the sphere, in contrast to (21) on the ball, does not correspond 
to the (^-function for any Laplacian. 

Finally, we note that (p d (s) does not have a pole at s = —1/2, which is im- 
portant during the construction of the heat-kernel coefficients and the functional 
determinants on M. It means that the log term is absent in the expansion (26), cf 
[1]. No further discussion of this term is given here. Some relevant comments can 
be found in the recent works by Bytsenko et al [19,20]. 

5. Residues and sphere heat-kernel coefficients. 

To apply (28) to the ball it is necessary to know the coefficients on the sphere. 
We use (27) on J\f with the appropriate (^-function, in this case either or £ . 

From (43), the residue of the modified coexact sphere (^-function at s = k/2, 
k E Z, is 

- Tk^>Y±- i y {d l 2j _ iKd-27-fc)! { p-J ) (48 » 
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where we have been able to use the results of our earlier work [1] for the residues of 0- 
form (-functions. The Dv are easily evaluated generalised Bernoulli polynomials. 

The corresponding heat-kernel coefficients are, for min ([n/2],p, d — p) = [n/2] 
with d — k = n, 

W'y / x . Mid+m) ^ (-iy^Dt?r 1} (d-2j-i\ 

\s d \ An / 2[p) [ft iJ r((fc + i)/ 2 ) ^ (d-2j-i)(n-2 3 y\ p-j )■ 

(49) 

A^ 2 (p) vanishes when n is odd and less than d because D„ is zero for odd v. 

A useful check is provided by the modified coexact coefficients A^ d ^ 2 = 
(—l) q Cp d (—q)/ql, q G Z + , for odd d. According to (45) these half-odd coefficients 
are non-zero, in agreement with the results of [21] and [22] for the de Rham coexact 
(-function. For this case, the only half-odd coefficient is the q = one, correspond- 
ing to C(0), and we know that the modified and de Rham coexact heat-kernels are 
related by a factor, 

K™ A (t) =e -((rf-l)/2-p) 2 r K deRham (r)5 ^ 

the expansion of which reproduces exactly (45). There are no half-odd coefficients 
for n < d, nor for p = when the zero mode is included. 

Similar statements hold for the 'total' coefficient, A n ^ . Because of the con- 
struction, (33), of Cp d+ 5 one should not, of course, expect a standard geometric 
form, except when p = 0. 

As a curiosity, constructing the 'total' sphere coefficient (c/(33)) by just adding 
equation (49) for p and p — 1, we have 

W J/ y + = (k _ ^((d+P/2) [ ^ ] 1V - 2^zfcr i) (d - 23 
\S d \ n ' 2 [ ' T((k + l)/2) l) (d-2j-l)(n-2j)\ \p-j 

(51) 

similar to the general form on the ball developed in the next section and also in 
agreement with the results of Giinther and Schimming [23] and Gilkey, e.g. [8]. 

Equation (50) can be used, in familiar fashion, to give general formulae for 
the standard, i.e. de Rham, coexact heat-kernel coefficients on the sphere, denoted 
£? n / 2 , somewhat simpler than those derived in [22] and [21]. Thus 

W w(p) M ~ N 



V ( 52 ) 
" w(p) M ~ N 



ns d / n _ A S d / >, 

%i/2^ - 2^ (M-NY N + 1 / 2[p) 
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where w(p) = ((d — l)/2 — p) and A^ 2 (p) is given by (49). For odd d, the only 
non-zero half-odd coefficient is Bfj 2 (p), as has already been remarked. 

Elizalde et al [21] give selected numerical values for the total B coefficients 
obtained from a general formula that involves a series of Hurwitz ^-functions and is 
slightly less convenient than the above. We have found agreement with their values. 

The connection between the two expressions is simply that the result of [21] 
follows from ours if the generalised Bernoulli polynomials are expanded in ordinary 
ones. From a numerical point of view this is unnecessary because there are efficient 
algorithms that give the generalised functions directly. 



6. General form of coefficients on the ball. 



We now come to the central part of this paper, the explicit construction of the 
heat-kernel coefficients on the ball. 

Firstly we can think of checking against the limited results of [2] for the coef- 
ficients for the form-valued Laplacian on manifolds with boundary. 

From Theorem 1.2 of [2] applied to the .D-ball, the total coefficient expressions 
for absolute boundary conditions are, (the + superscript is dropped in this section) 



(47T) d / 2 
(47T) d / 2 

\S d \ 
\S d \ 
\S d \ 



| CD) 



(?) 



1 



1 



2^(d+l) 



4>) 



A[ D \p) 



4/2 (p) 



i d 

d 



D 

P 

D 

p 

-6 



1 
1 

D — 1 
p-1 



+ 6 



384 



(13d + 2) 



288(d- 1) 



-4(29rf-26) 



D 
p- 



D-2 
p-2 

D - 1 
p-1 
D 

P 



(53) 



192(d- 1) 



+ 

-3 
-3 



Assuming, according to our general philosophy, that d is sufficiently large, it 
has been checked that the first three coefficients follow from our formulation using 
(28), (49) and various binomial identities including the recursion 



0- 



a-1 
b 



+ 



a-1 
b-1 



(54) 



The fourth coefficient will be confirmed later by a more powerful technique, and the 
series continued. 



12 



Although these results were derived on the basis of a restriction on p, no such 
restriction exists for the expressions (53) in [2]. The conclusion is that, once the 
coefficients are displayed as explicit functions of d and p, they can be continued out- 
side the restrictions, the reason being that their particular structure is guaranteed 
by the geometric formulation. 

These particular checks show that our general formulae are correct but it is not 
feasible to pursue the further evaluation of the coefficients in this way by algebraic 
rearrangement of the produced expressions, since it involves awkwardly repeated 
application of the binomial recursion and of other identities. A systematic approach 
to the evaluation of any coefficient is better provided by fitting unknowns in a 
general form, and this process will now be set in train. 

From basic theory (e.g. [8]) the geometric expression on a flat, bounded D- 
manifold, M., is, up to terms involving derivatives of the extrinsic curvature k, 

c{n){A^A° /2 {p) = I b n (D,p)J2(^(^)tr(^)...) (55) 

JdM n 

with 

c(n) = 2y / 7r, n even 
= 1, n odd 
= 2(rf+l)v^, n = 
and where n = (n 1 , n 2 , . . .) is a partition of n — 1. For convenience the n = term 
has been included although it is really a volume contribution. We seek the b n (D, p). 
For the .D-ball, (55) reduces to 

n-l 

c(n)(4n) d / 2 A° /2 (p) = \S d \ J2d ln K(D, P ) = \S d \ d k b£\D,p) (56) 

n k=l 

(n) 

where |n| is the number of components in the partition and b k , is the sum of those 
b n for which |n| = k, 

t\D,p) = J2K(D,p) 

n |n|=fc 

These are the only combinations that can be determined from working on the ball. 
The numerical multipliers b n (D,p) satisfy the binomial recursion 

b B (D,p) = b n (D -l,p) + b n (D - l,p- 1), (57) 

proved by crossing M. with a unit circle, [2]. This relation is what has become of 
the more familiar statement of dimension-independence for scalars. 
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Thus b n (D,p) can be expanded as a linear combination of binomial coefficients 
( p+6 ) f° r varying a and b. Since b n (D,p) vanishes for p outside the range to D, 
a must equal b and be nonpositive. Consequently the expansion reads 



n , 
b n (D lP )= J2 M n,J 

m=0 ^ 



D — m 
p — m 



(58) 



which is the boundary version of the Giinther and Schimming form [23] . The limits 
have been fixed by noting that they must be independent of both D and p and, 
therefore, can be set by considering the particular value D = n, for which there 
are n + 1 independent constants implying n + 1 terms in the sum. The trivial first 
coefficient, b = can also be used to pin down the limits to those shown. The 

conclusion is that, on the ball, the general form can be written, 



where Pm yd) is a polynomial of degree n — 1 in D. For n > 1 

n-l 

PL n) (d) = £ M n>m rfl n l = d k (60) 



n k=l 

where the M^ 1 ^ are constants. 

Our method is simply that, for fixed n, knowing the left-hand side for < p < n, 
(59) can be inverted for the Pm^ (<i), keeping d unspecified. The method employed in 
[2] uses the particular values b n {n,p) to invert (58). Our technique is more flexible 
in practice, although the analysis in [2] does bring out clearly the amusing fact that 
the general coefficient can be generated from just the values A^ 2 (p), < p < n. 

In detail, we firstly note that, for p = to n, the matrix of binomial coefficients 
on the right-hand side of (59) is triangular so that the inversion follows recursively 
by forward substitution, 



Pir ) M = £(!!^^) (m) _ m f ; > : ^P r(d) , m = „,...,„. (61) 



The driving coefficients A { n ^(m) (0 < m < n) will be determined from (28) 
as polynomials in d because, for given numerical values of p and n, the sphere 
coefficients, (49), are obviously such polynomials. 
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It should be stated that the ^-function on the ball is a standard one, and that 
the heat-kernel coefficients are geometric (leading to the general form (59)). In this 
case it is quite permissible to evaluate the polynomials in whatever parameter region 
is convenient and then continue, despite the fact that they involve coefficients on 
the sphere for which this might not be possible. 

Evaluation is a routine machine matter, the first outcome of which is (p) 

(n) 
m 



in (53). Some results are given below in the form of matrices of the constants M k 
in (60), 



2 13 _3 _1 
]V[( 3 ) = / 192 48 4 2 



13 

384 



29 
96 



(62) 



Mi 4 ) 



Ml 5 ) = 



V 



77 



15360 

235 
36864 

139 
61440 
2041 
737280 



4 


164 


16 


135 


" 315 


5 


1 


92 


74 


45 


105 


15 


1 


136 


26 


27 


315 


15 


77 




191 


960 




192 


263 




1475 


1440 




768 


1987 




1769 


15360 




1536 


3787 




347 



92160 



1536 



16 
3 



19 
6 

_47 

8 

157 
48 
__9_ 
16 



_ 15 

4 
55 
8 

_ 15 
4 



5 \ 
2 \ 

_ 11 

4 

3 
2 

A) 



(63) 



(64) 



A subscript has been added to indicate that these values are for absolute bound- 
ary conditions. Although any relative coefficient can be found from the absolute 
ones using duality, in order to obtain the general form the entire analysis should 
be repeated. A tactically better way combines these procedures and consists of 
using the recursive nature of the relative coefficients to write down a general form, 
precisely as in (59), together with the corresponding solution, (61), but now where 
the values of the coefficients for p = to n are determined by duality in terms of 
already evaluated absolute quantities. In any case, one rapidly finds that the first 
coefficient is the same, that the second is reversed in sign and that the remaining 
ones (up to n = 5) are contained in the matrices, 



Ml 3 ) = 



Ml 4 ) = 





5 


13 


3 




■(- 


192 


"48 


4 






7 


29 


3 


1 




384 


96 


4 


2 


8 


172 




16 


16 


189 


315 




5 


3 


11 


104 




74 




315 


105 




15 


8 


1 


116 




26 


8 


189 


315 




15 


3 



15 





109 


29 


193 


19 


15 


2 \ 

11 


15360 
1049 


320 
1247 


192 
1501 


6 
47 


4 
55 




184320 


5760 


768 


8 


8 


4 




47 


709 


1783 


157 


15 


3 




61440 


5120 


1536 


48 


4 


2 


V 


13 


2467 


325 


9 


5 


k 1 


147456 


92160 


1536 


16 


8 





(65) 



It is a matter of a few minutes by machine algebra to calculate larger matrices. 

Knowledge of these constants enables restrictions to be placed on the multipli- 
ers, b n (D,p), in the general geometric form of the coefficients, (55). 

Denoting the set of bj™\D,p), (1 < k < n — 1), by the vector b^ n \ the vector 
of binomial coefficients (0 < m < n), by and the set of powers d k 

(1 < k < n — 1) by d.( n \ equation (59) takes the matrix form 

^^HjQ(p) = dSOMWcO (66) 

and the restrictions are 

b W = M W c W. ( 67 ) 
For example, for n = 5 there are 4 partitions so 

1,(5) _ h -h +h /> (5) - h - h 

1 — 4 , 2 — 1 3 + 2 2, 3 — 1 2 )2 , 4 — l4 

and only these particular combinations are known using (67) and (64) or (65). It is 
clear that one can always unambiguously determine and 

We finally note that the n = 4 values at D = 4 agree with those in Table 2 of 
Moss and Poletti [24]. 



Generating functions. 

A generating function approach to the form order, p, could be introduced at 
this coefficient stage by defining 

A ( n %) = Ei-o-TA^ip) = '2 /2 ± {-ar{l - *) D - m P£Hd). (68) 

p=0 V A ) rn=0 

In particular the value a = — 1 is required (cf [2] Thm. 1.3), 

= ; w^TT^ E 2""PW(d) (69) 
c(n)r((d-l)/2)^ 

16 



which is a convenient way of organising and computing this quantity, at this point. 
There is topological interest in the case a = 1 and we see from (68) that 

= 0, n<D; ^(l| = ||gl)^M=^) = l 

(70) 

which provides a check of the Gauss-Bonnet theorem. If relative coefficients are 
used we find x(M,dM) = —1, which is correct. 



7. Value at zero. 



Special interest, both mathematical and physical, attaches itself to the par- 
ticular value C^ 4+ (0) which merits individual treatment, e.g. [24-26]. Further, as 
already said, these values play an important role in the algebraic computations of 
Blazik et al [2]. 

Generally, for absolute conditions, the coexact expression is 

C(o) = -^(-V2) + C^i(-i/2)) + \(C?(0) -C£i(0)) 

- E I (^"^ ( Res + C£i(*/2)) - ( - u a (p)) fc Res C^(*/2)) • 

(71) 

For Af = S d the residues are contained in (48) and the values of the ^-function 
found from (41) or (43) (there are no zero modes on the sphere in the sense that 
v > 0). 

An alternative evaluation consists of using the standard formula 

Ci 4+ (0) = A^ 2 (p)-Pi 4 (72) 

together with an explicitly computed general coefficient (59). The absolute Betti 
numbers of the ball are (3 Q = 1, the rest vanishing. 

The direct evaluation via (71) is the most efficient, numerically. Machine com- 
putation shows agreement between these two approaches and also with the values 
presented in [5], except that the p = numbers for d = 6 and 7 differ. 

It is a remarkable fact that in using (28), leading to (59) and thence to (72), one 
does not have to worry about the 'zero modes' in say (15), whilst such contributions 
are vital when evaluating directly from (71). 

Finally, using the inverse (22) and noting that Cd^O) vanishes, we reconfirm 
the topological formula X = ^ =0 (-l) p p£* = Ep=o(-^) P KM- 
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8. Form functional determinants. 



In this section we consider functional determinants for p-forms. As we have seen 
in the calculation of the heat-kernel coefficients, it will be enough to consider the 
determinant associated with the coexact (^-function (23), this being simply a com- 
bination of Robin and Dirichlet contributions. It is immediately appreciated that 
the results derived in [1], eq. (9.8) for Dirichlet boundary conditions and eq. (11.2) 
for Robin conditions, remain valid once the base (^-function there is replaced with 
the p- form base (^-function, eq. (14), with, in addition, 

and the related quantities, C^ +1 ( s ); { s i u )i if t ne zero modes are included as in 
(15). 

The coexact determinant for p > 1 combines Robin with Dirichlet and is de- 
termined, for absolute conditions, by the derivatives, 

C(0) =C^ +1 (0) + Cf_ + 1 1 (0) + (f(0,u a (p)) + In 2 (^(-1/2) +C^i(-V2) 

d d \ 

+ ResC^m) M^u.ip)) +2^ResC^ 1 (V2) D t (l) 

»=i i=l 

i odd 

d / i—1 

+ 2 Res C^(V2) ( M t (1, u a {p)) E V* 



i odd 



k=i 



f 1 M l (t,u a (p))-tM l (l,u a (p)) \ 



+ 2 E RfisC^(t/2)(M 4 (l,« a (p))X;V* 

«=1 ^ fc = l 

i even 



i-l 



L 



+ / dt 



t(l-t*) J 



+ 2±Res C^/2) (^(1) £ l//c + £ dt ^^f ' 

(74) 

The M i and Dj are known polynomials associated with the asymptotic expansion of 
Bessel functions and, to avoid repetition, we refer to [1], and to references therein, 
for more details. Some information can be found incidentally in Appendix D. 
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For p = a small extra consideration is necessary. As we have mentioned, 
absolute 0- forms are Neumann scalars and the Robin parameter is then u a (0) = 
— (d — l)/2. Looking at (^f (s,u) in eq. (73) it is seen that, for u = —v{p), one 
encounters a branch cut occasioned by (incorrectly) including the true zero mode 
for Neumann conditions. The technical reason is that the asymptotic expansion for 
these specific Robin parameters is slightly different from the others. 

The easiest way of taking this into account is to subtract the contribution in 
(73) for u 7^ —v(p), then take the limit as u — > —v(p) and finally to add the correct 
contribution for u = —v{p). 

For Neumann conditions the end result is 

C'^(O) =C^ +1 (0) + ln2fc;f(-l/2) + 2 £ Res ${i/2) M,(l, « o (0)) 



i = l 

i odd 



i-1 



+ 2 Res^/2) (M t (l,u a (0))J2 1 / k 



i=l X fc=l 

i odd 



M l (t,u a (0))-tM l (l,u a (0)) \ 

1 ^ " t(i-t*) ^ ; (75) 



+ 2 ResCo^/2) M { (l, Ufl (0))^l/fc 



<=i x fc=i 

i even 



+ / 



M t (t 7 u a (0))-t 2 M z (l,u a (0)) \ 
t{l-t*) J 

+ lim fCo V (0,«)+ln((d-l)/2 + «)) +ln(d+l). 

u— ► — (<i— 1)/2 V / 

The relative, Dirichlet expression is given in [1]. 

Equations (74) and (75) are expressions on the generalised cone and we apply 
them now to the ball, for which all quantities have already been treated, apart from 
C^ +1 (s) (C^ +1 ( s )) an d ( s i u ) (Cp ( s , u ))- For these remaining ^-functions one 
immediately finds, along the lines of [1], 

t P V+1 (*)= E ( m ~ 1 )(B(sAd + ^/2\l m+1 )+C R (s;(d + 3)/2))6 pd 

m=p+l VP/ 
d 



+ E ( d ™ ^ 1 )c B (a»(d + 3)/2|l m+1 )+C Jl (a;(d+l)/2))(5 1J o 

m=d—p 

(76) 
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and 



Ct r (s,u a )= f: ( m v 'W^+ilU 

m=p+l VP/ 

+ ± (r^y^^- (77) 

m=d—p 

+ d p0 {(d-l)/2 + u)- s + d pd (d+ir s - 

The contribution of the zero modes is visible and (^ +1 (s), (s,u a ) need not be 
stated. It is clearly seen that the limit u — > —{d — l)/2 in (75) is well defined 
because the logarithm is cancelled by the last term in (77). 

Using (76) and (77) in (74) and (75), the determinants emerge as derivatives 
of the Barnes ^-function at s = 0. These may be expanded using Stirling numbers, 



C B (s,a\l m ) = jr( m + ™ X Va + n) s 

n=0 V / 
™ ( — l^k+m 

= g (Jl) l ( m - fc )! B a(°)C B (. + l-*.a 



to give derivatives of Hurwitz or Riemann ^-functions, [27] p. 433. The derivatives 
presented in Appendix B were obtained in this way. It is easy to find the expression 
in any dimension d and for any value of p. 

The presentation of our results is a little more explicit than the corresponding 
ones of [5] in that only derivatives of Riemann C-functions are involved, there being 
no integrals over r or $ functions. In fact, the integrals in [5] can be done, bringing 
exact agreement between the two sets of expressions. This agreement is elaborated, 
in its essentials, in Appendix C. 

It should be noted that, as a consequence of an earlier error in [9] for Neumann 
0-forms, a log 2 term should be added to the results in [5] for p = and d > 2. 



9. Spectral properties of spinors. 

For the sake of completeness, the method of [1] is now employed to put the 
results of the earlier work on spinors, [28], on a secure and systematic footing. For 
any undefined quantities and for further explanations, refer to this paper. 

Spectral boundary conditions. 
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In [28] it was shown that spectral boundary conditions lead to 



J n+D/2 _ 1 (/c) = 0, neZ, (79) 

with the degeneracy 

d(n) = d s ( n + D - 2 ). (80) 



n 

Having in mind the treatment of the scalar field in [1] we define the spinor base 
^-function appropriate for the ball, 

oo 

C^(s) = J2d(n)(n + D/2-l)- 2s . (81) 

n=0 

This choice guarantees that all formulas in terms of ^ m section 3 of [1] remain 
valid. Also the polynomials D^t) defined there are unchanged. In particular, 
eq. (4.8) is still true, the only difference is in the base (-function used for the 
calculation of A^, , the base heat-kernel coefficients. 

n/2' 

Let us continue with the treatment of (_^ (s) for the ball. It has the especially 
simple form 

(js(s) = d s ( B (2s,(d-l)/2\l d ). (82) 

To obtain the heat-kernel coefficients using eq. (4.8) of [1], we need only the residues 
in (82). Using the known values of the Barnes (-function in eq. (4.8) of [1] one 
regains all the results of [28] for the heat-kernel coefficients including the conjecture 
for their general form. 

The determinants come out similarly. Due to the identical structure of the 
present problem, eq. (9.8) of [1] remains valid once the base zeta function there is 
replaced by (81). The only 'unknown' quantity is 

CAr+i(s) = d s C B (s,(d+l)/2\ l d+1 ). (83) 

Using (78) this can be expanded in Hurwitz-Riemann (-functions, and complete 
agreement with the results exhibited in [28] is again quickly and economically found. 

In Appendix D we prove the curious identity for the spectral heat-kernel coef- 
ficients, 

4d-2)/2( D ) = °> (84) 

observed numerically in I, the significance of which has not yet been discovered. 
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Local boundary conditions. 



For local boundary conditions a few more technical details are necessary. The 
eigenvalue condition is, 



Jn+D/2(k)-Jn+D/2-i(k) = 0, neZ (85) 
with the mode degeneracy 

which is half that of the spectral case. Defining the base ^-function as in (81) we 
have 

Cx{s) = \d s C B (2 8 ,D/2-l\l d ). (87) 
The polynomials describing the asymptotics are a bit different, however, and read 

A(*) = E*i,a* i+ °, 
a=0 

which leads to a slightly different form of the A^s). These are found to be (one only 
has to refer to [29] eq. (2.18), and replace the ^-functions by the base ^-function), 

A 1 (a)= V r( !~ 1/2) CAr(s-V2), 
~ 1V ; 2v^F r(s + l) wv 1 ' 

MS)= ~2^ T(s + 1) C " (S) > ( 88 ) 

MS) ~ T (s)^ {S + l/2) ^ a T((z + a)/2) • 
Due to the different A (s) and A^s) the heat-kernel coefficients this time read 

= ri %D-T ) Baa ^ ((D - n - l)/2) 

v- 1 ^ > r((D-?i + z + a)/2) 

i=l a=0 ~~ ^ 

It would be possible to rewrite the right-hand side in terms of the spinor heat- 
kernel coefficients on the base, leading to an equation like (28). 
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For the ball, using (48) for the residues of (j^(s), (87), all results of [28] are 
immediately reproduced in an organised and rapid fashion. Let us therefore proceed 
to the determinant. 

As a result of (88) one has 

^(0) =2 (In 2 - 1) Cat (-1/2) - ("1/2) , 

1 



4,(0) =^2^(0) --c^(o), 



2i 



A'm = -J2 *i,a [PP Cm (i/2) + 7Res Cm (*/2) + 



(90) 



a=0 



Res 0v(i/2)^((a + i)/2) , 



with characteristic differences to the spectral case. Further contributions come from 

oo 

Z'(0) = ^2d(n)Z v, (0), u = n + D/2-l, (91) 

and 



n=0 



D-l 



Z u '(0) = 21nr(z/ + l) + 2z/-2z/lnz/-ln(27rz/)+ ~ ^ 



71=1 



(92) 



This follows from the results preceding eq. (3.1) in [29]. 

Formally (92) is twice the scalar case. It also shows that D n (l) = 2( R (—n)/n 
and so the contribution Z'(0) is known from the scalar field calculation. Adding 
up, everything fits nicely and the final answer reads 



C| o /2 '(0) = ln2(2Csr (-1/2) + Cat (0) ) + 2& +1 (0) 



D-l 



+ 41n2^^-^Res (*/2) 



i=l 



(93) 



D-l 



+ 2 ]T Res Ca^/2) 



i=i 



2C fl (-Q gl , /"^ D t (t) -tD.il) 



k=l 



t(l-t 2 



Specialising to the ball, all results of [28] are thereby established again in a system- 
atic and general way. 
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10. Conclusion. 

Our essential results are, firstly, the expression (28) for the coexact heat-kernel 
coefficients on the generalised cone in terms of those on its base. Secondly, (59) gives 
any ball heat-kernel coefficient, with explicitly computable polynomials, (60), and 
leads to the restriction (67) on the general coefficients. Thirdly, we draw attention 
to the formula (74) for the functional determinant and finally to the expression 
(43) for the modified p-form coexact (function on the sphere as a finite sum of 
(conformal) 0-form functions. 

To save excessive analysis, we have concentrated on the simplest generalised 
cone, i.e. the ball. As a continuation, more complicated bases, A/", can be consid- 
ered. In the earlier paper, [1], we looked at the case when the base was a sphere of 
non-unit radius, say a, calling this a monopole. Then, in place of (35) for the order 
of the Bessel functions, we have 

A^P) = ^(l + (d-l)/2) 2 - ((d- l)/2 -p) 2 (± - 1) (94) 

which is not a perfect square unless p = (d — l)/2. We therefore expect that a pole 
at s = in (^(s) will arise even for the case of the de Rham Laplacian. However, 
unlike the scalar (p = 0) case, it does not seem possible to achieve a perfect square, 
and hence to avoid the pole, simply by attaching a term to the de Rham Laplacian. 
This obstruction is probably related to the fact that there is no direct way of making 
the equation conformally invariant except for the special value p = (d — l)/2. It 
is possible to use Branson's operator [30] but we would then lose contact with the 
body of work associated with the de Rham Laplacian. 

The lack of conformal invariance can be appreciated in the structure of the 
coexact (^-function on the sphere, (43) , which is a combination of improved 0-form 
(^-functions each of which is conformal in a different dimension. 

Further possible choices for Af include hyperbolic spaces and factored spheres 
for which there would be more topological and analytical excitement involving, as 
they do, multiple connectedness and conical singularities. 

Regarding spectral spinors, our previous method involved fitting a conjectured 
general form for the coefficients using specific numerical values. The present tech- 
nique establishes both this general form and the actual expressions in one operation. 
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Appendix A. 

In this appendix we derive the form (43) for the coexact sphere (^-function. 
The most important property for the proof is the 
Lemma. 

]T ( m ~ l ) C B (2s; (o+ l)/2 | 1J = - fl ( m ~l) (a - l)/2 | lj 

m=p+l \ " / m=p ^ ' 

d-1 



+ [ p JCs(2a;(a-l)/2|l d ). 

(95) 

Proof. 

Rewrite the left-hand side of (95) using the difference relation, [31], 

( B (s; (a + l)/2 | lj - ( B (s; (a - l)/2 | lj = (a - l)/2 | l m _J (96) 

and then combine terms using the binomial recursion (54). 

Now we are set to show the following recursive property of the p-form ^-function 

(41) 

Lemma. 

cfw+rrw= ( <i ; 1 )co sd w+^, 1 (^)" 2 '^,o(^)" 2 ' m 

for <p < (d- l)/2. 
Proof. 

Starting from (41) and (15), use of (95) gives first 

Cf («) = - E (Tli) <^ 2s; ( rf " I ^) + ~ *) (d - i)/2 | i d ) 

+ e (^^Jc^^+d/zhj+v^)" 2 *- 

m=d—p ^ ' ^ ' 
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Separate the m = d term from the last sum, then 



Cf(s) = ( d (Ch(2*; (d+ l)/2 | l d ) + Cs(2 S ; (d - l)/2 | l d 
" E (Tli ) Cb(2«; (d - i)/2 | i m ) + 5 P , (^) 

+ E ( d T~^i)c B (2«;(d+i)/2|iJ. 

m=d—p 

Next, combine the m = d — 1 terms of both sums using (96), 

<f («) = ("l *) Cf M - E ("I*) C»(2.;(<l- D/2 I i,J 
+ E (^"_ 1 1 )cB(2 S ;(d+i)/2|iJ 

m=d—p 

- : i ) cb(2*; (d - 1)/2 1 i d _ 2 ) + 5 P)0 (^) ~ 2S . 

Apply (95) once more in the second sum and then the binomial recursion, (54), to 
reach Lemma (97). Finally, applying this Lemma p + 1 times, for p < (d — l)/2, 
produces the important relation (43) in the main text. No doubt this also follows 
from a direct analysis of the degeneracies. 

A comment on these equations is perhaps necessary. Setting p = in (97) we 
must recognise that C-i ( s ) equals ((d — l)/2) s for consistency. This follows 
easily from (41). In general, C-p( s ) = ((d — 1 + 2p)/2) 2s . 



Appendix B. 

In this appendix are listed some selected derivatives at zero of the coexact zeta 
function for absolute boundary conditions and with M the ball. 
In d = 2, 

C - = _g_^ + log3 _^ + ^l) +m 

C- = -l + ^-^) + 3C'(-i)-C'(0) 

/M _ 3 log2 3C ; (-2) C(-l) 
^ 2 32 12 4 2 
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In d = 3, 

M _ 1213 151 log 2 CM) C(~2) , 13CX-1) , , Vn , 
Co " 4320 + 90 + 3 + 2 + 6 + C(Uj 

>,_ 5989 19 log 2 C(~2) 3C'(-1) 

Cl "10080 " ^O - + C { ~ 3) + ~~2 2 C {U) 

^=-^ + ^r + c(-3)-^-^ +2 c(o) 

^ = J73_ log2 C(-3) _ CV2) C ; (-l) 
^ 3 30240 90 3 2 6' 

In d = 4, 

, M= 25381 17 log 2 5C ; (~4) 23C'(-3) 

^° 46080 2880 g 64 48 

47 C '(-2) 108f(-l) 

32 48 S W 

M = J3803_ 77 log 2 _ _ 5C(-4) 19C'(-3) 

^ 11520 720 g 16 12 

17 C -(-2) _ 25C'(-1) 

8 12 s v ; 

M = ^09_ _ 863 log 2 _ 15C'(-4) 15C'(-3) 

^ 2 2560 480 32 8 

_3C( z 2)_21C-(-l) 

16 8 S V ; 

M = _ _2509_ 77 log 2 _ 5C(-4) llC'(-3) 

^ 3 11520 720 16 12 

-^ + i^)-C'(0) 

M = _4J_ 17 log 2 _ 5C(-4) 7C(-3) _ C'(~2) _ ('(-!) 
U 9216 2880 64 48 32 48 

The (^-functions on the right are Riemann ^-functions. 

As mentioned in the text, it is possible to obtain the determinants in any 
dimension d and for any value of p without difficulty. 



Appendix C. 

In this appendix we outline a method for the calculation of integrals of the type 

,6 

/ dww^logTiw), (98) 

J a 

needed in order to see explicitly the agreement of our results with those of [5]. 
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We start from the usual relation 

^C(8,w) = -8<Z(s + l,w), (99) 

valid for any modified (^-function, such as (^f (s) or the Barnes (^-function, but here 
applied to the Hurwitz ^-function. 

The differentiations are continued to give, 

((s,w) = (-l) n s(s + l)...(s + n-l)((s + n,w) = f n (s)((s + n,w). (100) 



Now differentiate with respect to s 

Qn 



and set s = —n 

Qn 



C'(s, w) = fMC(s + n,w) + f n (s)('(s + n, w) (101) 



C(-n, w) = /;(-n)C(0, w) + f n (-n)C(0, w) 



= /;(-n)C(0, u;) + / n (-n)(lo g r» - log v^F). 
Next, multiply by and integrate 

/•b an 
,n-l u»rL^ _ / J„..„.,n-1 ^ a// 



/ n (-n) y 1 log !?(«;) = y dw«; n £'(-71,11;) 

+ ^t^(b n -a n ) log V 7 ^- fn(-n) [ dww^C&w). 

(103) 

Define 

/b an 

and integrate by parts to get the recursion 



n-l & 



/ ( s ) = w n ~ 1 -^ r('(s,W 



(n -!)/„_! (a) (104) 



which is trivially continued using ^(s) = C'(s, b) — £'(s, a ) 
From (101) we have 



am 

w m —('(-n, w) = w m (/^(-n)C(-n + m, «;) + / m (- n )C'(-n + m, w)) (105) 



<9u; m 
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where 



m 1 

U-n)=ml( n ), fL(-n)=ml( n )j2 

\m/ Vm/ ^-^ m — r 

k=i 



n — k 



(106) 



Using the recursion (104) n times, and also (106), the following result is found for 
the integral (98), 



(107) 



/ (1 1 I 1 

dw u,"" 1 log T(w) = - log( 27 r) + T w n ~ Z—~^ n+1 E Z 
^ L k=i J fe=i 



i=i 



(-1) 



z-i 



n — 1 
/-l 



n—l 



c(-i,w)-a-i,w)j2 



k=i 



k+i 



The simplest case, when n = 1, is given in [32] and an integral related to (107) for 
n = 2 can be found in Nash and O'Connor, [33] eqn.(C.30). 

If one of the limits in (98) is zero, the last sum in (107) has to be interpreted 
as the limit as w — > and, in this case, we use 



C(s,w) = w 3 + C(s,l +w) 



(108) 



to see that 



lim w n 1 



n—l 1 



k=i 



k + l 



= 



for / < n, whereas 



lim('(-n, W )=C'(-n) 



for n > 1. 

Let us now give some examples of the general formula (107). To compare with 
Co(0) in d = 3 of [5], we need 

<iu; (w — 1) logT(w), 

which is obtained by setting n = 1 and n = 2. For n = 1, we get the well known 
formula ([34], 1.9.1 (18)) 



rf«;logrH = ilog(27r) = -C / (0) 



(109) 



which is the essential part of Raabe's formula. 
For n = 2 



/ 

Jo 



dww\ogT(w) = ~- ^C'(O) + C'(-l)- 



(110) 
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The integral term in [5] is therefore 

r 1 1 

2 J dw(w-l)\ogT(w) = -- + 2C / (-l)+C'(0). 

In other dimensions, generally a and b take integer or half-integer values. Using 
(108) the integrals can always be expressed in terms of derivatives of the Riemann 
(■-function. We meet some known integrals such as( [35,32]) 

/3/2 3 13 ^ 1 

dw logI» = -- - - log 2 - -C'(-l) - -C'(0) (HI) 

and some lesser known ones, e.g. 

J* /2 dwwlogr(w) = -\ - ^log2 + |c'(-2) - ^C'(-l) - §C(0) (112) 

generalisations of which can be found using (107). 

Finally we point out that a generalisation of Raabe's formula to the multiple 
T-function is derived by Barnes [31] using the same ingredients as in the above. 



Appendix D. 

The identity (84) will now be proved. In terms of the spinor base (-function 
it reads 

d-2 

Res Cat (V2) -Res Cat (1) - 2 ^ Res (1 + i/2) D[{\) = 0. 

i=l 

Using knowledge of the residues in terms of the Bernoulli functions, [31], one has 
flgi((J-l)/2) Bff 2 ((d-l)/2) 

(d-iy. (d-2)! 

d-2 „(d) „, y LL6 ) 



+ 22^ 1J ^UJ (1 + i)!(d _2_0!' 



, ,.(d-2-i)! 

which is a standard recursion formula for the generalized Bernoulli functions as will 
now be shown. 

We first note that -D^(l) = (— 1)*/2. This value is obtained by comparing the 
small z approximation of the Bessel function I v {yz) with Olver's asymptotic form. 
Thus, working always to order z 2 , 

In />*) ~ v]n (uz/2) - \nT(l + u ) + ]n(l + ^-^{yz) 2 ^ . 
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Expanding the log, 



In ( 1 + —±—(vz) 2 ) ~ ^ 2 ^r^ = ±z 2 i/ 
V 4(z/+l) v J 4 1 + l/z/ 4 ^ v 



and also 

lnlYl + z/) ~ i/lni/ - z/ + - \n2nv - V ^IZ^ -fc. 



k=i 

The Bessel asymptotics are 



11 / z \ °° 

\*I v {yz)~-- \u2tvu-- ln(l + z 2 ) + i/^>/l + * 2 +ln ( ^-|— )j+^ 



n=l 



~ -- In 2™ - -z 2 + + -z 2 + Hz/2)) + ^ n 

n=l 

Comparing, we obtain the quoted value of D' n (l), and also that of D n (l) met 
in our previous work. The higher derivatives can clearly be found by carrying the 
expansion further. 

Substituting these values of D' n (l) into (113) yields the condition 

E( d 7 1 )^ 1 _ i ((d-l)/2)) = (114) 
which is a special case of the recursion relation ([36] eqn.(ll) pl61) 

j^( y )^B^} J (x\u)=uvBt\ 1 \x\u'), « = (115) 
j=o ^ J ' 

involving the general Bernoulli function. Setting u> = 1, v = d — 1 and n = d, this 

becomes 

d-i 



E ( d ; ')*-;(') = (''- Iptf-Vto 
i=o v J y 



(116) 

= (d-l)(x-l)(x-2)...(x-d + 2). 

If (i is odd, (<i — l)/2 is integral and, since 1 < (d — l)/2 < d — 2 for d > 1, the 
right-hand side of (116) vanishes for all odd d when £ = (d — l)/2, as required to 
show (114). 

For even d, = 2q, (116) allows one to give the general form of the coefficient 

^-2 )/ 2^) = {-iy{ q -\){{ci-\)...^ , £> = 2 ? +l. (117) 
A similar result should hold for the monopole. Put to = a in (115). 
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